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Abstract 

> : 

I The thermal and statistical properties of hadronic matter under some ex- 

treme conditions are investigated using an exactly solvable canonical ensemble 
I model. A unified model describing both the fragmentation of nuclei and the 

thermal properties of hadronic matter is developed. Simple expressions are 
^ ' obtained for quantities such as the hadronic equation of state, specific heat, 

c ■ 

compressibility, entropy, and excitation energy as a function of temperature 



^ I and density. These expressions encompass the fermionic aspect of nucleons. 



such as degeneracy pressure and Fermi energy at low temperatures and the 
ideal gas laws at high temperatures and low density. Expressions are devel- 
oped which connect these two extremes with behavior that resembles an ideal 
Bose gas with its associated Bose condensation. In the thermodynamic limit, 
an infinite cluster exists below a certain critical condition in a manner similar 
to the sudden appearance of the infinite cluster in percolation theory. The 
importance of multiplicity fluctuations is discussed and some recent data from 
the EOS collaboration on critical point behavior of nuclei can be accounted 
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for using simple expressions obtained from the model. 
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I. INTRODUCTION 



One of the goals of nuclear physics is to characterize the thermal properties of nuclear 
matter, described by its equation of state. Investigating heavy ion collisions offers one 
experimental method of probing nuclei away from the typical low temperature and density, 
and a number of experiments are underway studying the breakup of nuclei, contributing 
greatly to our understanding Several theoretical investigations of these properties have 

been instigated, thermal |]7|-pi| , statistical |[n|-|T3| transport [T^-|5D[ and percolative [^1|HH 
models all offering insights into the properties of fragmenting nuclei. 

In a previous set of papers |p7|-p0| an exactly solvable canonical ensemble model was 
developed for studying the multifragmentation of nuclei induced by high energy collisions. 
Detailed properties of the cluster distributions such as inclusive yields and correlations were 
explored. In this paper, we study the statistical and thermal properties of nuclei using the 
same model, extending results presented in [^]. This model gives a unified description of 
multifragmentation phenomena and the thermodynamic properties of hadronic matter, such 
as its equation of state and the nuclear compressibility. Simple expressions for these quanti- 
ties are obtained which extrapolate or connect between a low temperature nearly degenerate 
Fermi gas to a high temperature ideal Maxwell-Boltzmann gas of nucleons. Importantly, 
the effects of the fragmentation degrees of freedom on the thermodynamics is manifestly 
included. Questions related to very high temperatures and particle and anti-particle pro- 
duction will be developed in a subsequent paper. 



II. MODEL STUDIES 

Given the partition function Z^(V", T) for a system of A nucleons at a given volume 
V and temperature T we can determine the thermodynamic functions by taking various 
partial derivatives, i.e. the internal energy and pressure are given 

P = ksT {-^XuZa) ■ Additional partial derivatives lead to such important quantities 
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as the specific heat and isothermal compressibihty, Cy = = ~ iXw/T^ 

thermodynamic potentials such as the entropy and Gibbs free energy can be related to these 
quantities, S = U/T + kBlmZ^, G = —ksT In Z a + PV . Our interest of course is not in 
the partition function itself, but rather in these derived functions, and as such the partition 
function is merely a means to that end. 

This section constructs a partition function for fragmenting nuclear matter. It does so 
by investigating the thermodynamic functions at both low and high temperatures and densi- 
ties. At low temperature/high density, the nucleons coalesce into a single nucleus which can 
be modeled as a degenerate Fermi gas with Skyrme interactions. At high temperature/low 
density, the nucleons evaporate into a Maxwell-Boltzmann gas of individual nucleons. With 
these results firmly in mind, a statistical fragmentation model is introduced whose parame- 
ters depend on the thermodynamics in a way consistent with the high and low temperature 
limits, yet allowing more complicated fragmentation patterns at intermediate temperatures. 
RHS always diverges in the infinite limit The remaining freedom in the model is elimi- 
nated by assuming that the cluster yields follow a power law. Interestingly these minimal 
assumptions lead to a model surprisingly rich in physics as we shall see in the next section. 



A. Low temperature behavior and ideal Fermi gas laws 

At zero temperature, the nucleus can be treated as a degenerate Fermi gas, with Fermi 
momentum pp given by 

gsj^^pl = h^P , (1) 

where gs,i = 4 is the spin-isospin degeneracy factor and p = A/V is the density. At normal 
density po = A/Vq = 0.17 nucleons fm~^, nucleons are non-relativistic and the Fermi energy 
and momentum are related by ep = p\/2m. From the above expression for the Fermi 
momentum, we see that ep oc p^/^. At low temperatures, the total kinetic energy for such a 
Fermi gas is given by 
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However, the nucleons do interact. Interactions between the nucleons can be taken into 
account by using a density dependent Skyrme interaction |33,^|. In this approach, the 
internal energy of a nucleus is given by 

Ui,T = EK-A(^ao-^^-a3(^£j j , (3) 

where oq and are Skyrme parameters (for simplicity, temperature independent) deter- 
mined by fixing ay(p) = —UuriT = 0)/A the binding energy per particle at zero tempera- 
ture 

«v(p) = -?eF(p) + oo ( — ) - as ( — ) (4) 
5 \PoJ \poJ 

to the empirical value for ay at p = po (about 8.0 MeV/nucleon at zero temperature) and 

having this as the maximum of ay- The density dependent repulsive term should appear 

with a higher power of p than the attractive part, so that the nucleus does not collapse, and 

thus requires a > 0. Commonly used values for a are a = 1 (three body interaction) and 

(7 = 2/3 (finite range term) ||35|| . 

Having specified the internal energy, we can now determine any other thermodynamic 

function. The specific heat Cy = (f^)^? and entropy S = Cy/T' dT' are therefore 

=(!-] =A-(^] (5) 
Ub/lt Ub/lt 2 [ep J ' ^ ^ 

The partition function can now be determined from Za = expK^T — U)/ ksT}, which gives 

ZA = y^, (6) 

where 

f ay , ksT) 

and £0 = ':^^F- /,From this partition function it is easy to determine the pressure and 
compressibility. 



2 _ Stt^ fkBTV 



Plt = -P£f{p) 1 + 



5" \ 12 \ ep J 

-a^p-^ + (1 + u)a^p - , (8) 
1 \ 2 . TT^ fksTV^ 



(-) -i'-^ 1+ 



12 V £f 



2 / \ 1+<T 

-2ao- + (l + (7)(2 + (7)a3pN^ . (9) 
Po \PoJ 



B. High temperature, low density ideal Maxwell-Boltzmann gas laws 

At high temperature and/or low density, but below meson and particle- antiparticle pro- 
duction thresholds, the hadronic properties are that of an ideal Maxwell-Boltzmann gas, 
with 

Uht = lAksT , (10) 

"--B / HT ^ 

Pht = pksT , (12) 

— ) = pkeT , (13) 
Kr/ HT 

Pt =Ai.y^^gsX, (14) 



where the entropy S is given by the Sackur- Tetrode law and At = /i/\/27rm/csT. We expect 
that nuclei vaporize into individual nucleons at ksT ^ oy, the binding energy per particle. 
The partition function consistent with these thermodynamic functions is Za — x^/A\, where 

X = y/A^. 

C. Hadronic matter at moderate temperatures and densities and ideal Rose gas- like 

structure 



In the region between the ideal Fermi gas limit and the ideal Maxwell-Boltzmann gas 
limit, the effects of the fragmentation of the initial nucleus must be taken into account 
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when considering thermodynamic issues. To account for these effects, we employ a model 



developed initially to study multifragmentation phenomena . This model has the correct 
high and low temperature limits, as determined in the previous sections, but has features 
similar to that of a Bose gas in the intermediate range. 

To describe the fragmentation of a nucleus into all possible modes of breakup, a weight 
is given to each possibility. The weight chosen is 

where x and y are functions of the thermodynamic variables (V,T), (3k is the cluster size 
dependence of the weight and n = (ni,n2,...) is the fragmentation vector, with the 
number of fragments with k nucleons such that J2k kn^ = A. With this choice of weight the 
canonical partition function Za = J2n W{n) is a polynomial in x, y, given by 

A 



ZA{x,y)=Y.zP0)x-y^-^, (16) 

m=l 



with m = J2j the multiplicity. Given the partition function, all thermodynamic properties 
of the model can be obtained, as well as ensemble averages. For example, the mean number 
of clusters of size k is 

xy^-^ ZA-k{x,y) . 

{rik) = —5 ^7 r ■ (1^) 

f3k ZA{x,y) 

The partition functions themselves can be obtained using a recursive procedure defined by 
the constraint J2k k{nk) = A, which gives 

1 ^ xy''~^ 

Za{x, y) = -rY. k——ZA-k{x, y) , (18) 
A (3k 

where Zo{x,y) = 1. The whole procedure is easily implemented by computer. 

The parameters x, y determine the thermodynamic aspects of the models and as such 
need to be determined. They are most easily determined by considering the partition func- 
tion at high and low multiplicity. When the multiplicity is large (m) ^ A, and Za ~ Za^'^x'^. 
Clearly this is the high temperature limit, and x is simply the ideal gas x introduced in sec- 
tion [II B| , given by 



^ = ^ (19) 

and involves the volume of the system V , the thermal wavelength A-r, and the dimensionality 
of the system d. The term in the weight arises from the thermal motion of each fragment. 
Since the overall motion is zero, this should be replaced by to reflect the center of 

momentum constraint. When the multiplicity is small (m) ^ 1, and ~ Z^^xy"^^"^. 
Assuming that the x parameter is removed by the conservation of momentum constraint 
describe above, this shows y to be given by Eq. (^. A correction to the internal excitations 
however needs to be made at high temperatures to reflect their finite lifetimes. Koonin and 
Randrup ^\ argue that a simple cutoff is effective, which results in 

Here ay is the binding energy per particle, Eq is the level spacing parameter (related to the 
Fermi energy) and Tq is the cutoff temperature. So the term y^""* in the weight arises from 
the binding and internal excitations of each fragment. 

The only parameter now to determine is jSk- Since the mass yield distributions are often 



well represented by a power law |3^, we choose Pk simply to reproduce this important 
experimental fact. If j3k = k'^ , then in the grand canonical limit (uk) ~ k~'^. Typically r is 
found between two and three, and we have chosen r = 2.5 as representative for calculations 
made in this paper. 



III. THERMODYNAMIC PROPERTIES OF HADRONIC MATTER 

In this section we consider the various fundamental thermodynamic functions of frag- 
menting nuclei, namely the pressure, compressibility and Gibbs free energy (section |III A| ), 
the internal energy and specific heat (section [III B[ ), and the entropy (section |III (J| ). Inter- 
esting features of the functions are explored, although a detailed examination of their origin 
is postponed till the next section. A brief accounting of such properties has already been 
given [ig. 
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Since the partition function can be expressed succinctly as a function of x, deriva- 
tives of In Zj^ with respect to these variables will appear and reappear in computations 
of thermodynamic functions. Fortunately for this model, derivatives with respect to x, y 
yield expectation values of the multiplicity, i.e. x^lnZ^ = (m), y^lnZA = A — (m), 
(^x-^^ IuZa = (y^) InZ^ = Var(m) = (m^) — (m)^. A general derivation of these identi- 
ties is included in appendix |^. 

A. The hadronic equation of state and a van der Waals like structure 

The hadronic equation of state is of much interest since the pressure and incompress- 
ibility reflect the behavior of matter at fixed temperatures and varying volumes. A simple 
expression can be obtained for this equation of state by using the P = k^T^ In Z^ and the 
partition function of Eq. (|16D as determined from the weight in Eq. (|15D , giving 

"-('"^PHX+fl-MUx, (21) 



where (m) = is the mean multiplicity, Plt, -Pht are the low and high temperature 

limits defined in sections [11 A| , [11 B| and = Plt — fp^^o (^^) t+To simply the cutoff 



corrected low temperature pressure. The above equation shows the importance of the mean 
multiplicity in the behavior of P. For low x/y, (m) ^ 1, and the pressure reduces to 
the Fermi gas and Skyrme pressure Plt- At high x/y, (m) ^ A and the ideal Maxwell- 
Boltzmann gas pressure Pht in P dominates. Thus the result of Eq. (pl[) connects the two 
extremes in a simple analytic way. Figure |l](a) plots the pressure for several temperatures 
for (3k = k^'"^- 

The low temperature component of the pressure is quite interesting. It is responsible for 
the pressure rising as the volume increases over a range of volumes at low enough temper- 
atures, a feature characteristic of a van der Waals gas or a liquid-gas phase instability [§. 
It arises due to the nuclear interactions. To see this phase instability more clearly, we can 
calculate the isothermal compressibility from the equation of state, since = (p§^)t, a 
little effort gives us 



1 ("^) T-> , f ^ ("^)\ "^-^LT 



Pht + ( 1 - ^ ) P- 

(22) 



Kt A \ A J dp 

Var(m) (Pht - PIt? 



A Pht 

At large values of x/y, (m) A, Var(m) and ^ ^ Pht? the ideal gas limit. At low 
x/y, (m) ^ 1, Var(m) ^ and PLti ideal Fermi gas with Skyrme interaction limit. 

In between, there is an effect due to the fluctuation in the mean number of clusters. The 
behavior of the compressibility is shown in Fig. |l](b). From the figure, the phase instability 
appears first at p = O.Spo? T = 24 MeV, the point where and its derivative are both zero. 

Another method of seeing this instability is to consider the Gibbs free energy G = 
—ksT In Za + PV . Figure ^ shows that G vs. P is multivalued below a critical temperature, 
a clear indication of a phase instability. 

In nuclear physics, one usually refers to the incompressibility k = which if ep ~ 37 
MeV gives at T = 

K = 96 + 144aMeV (23) 

For (7 = 1, this model k = 240 MeV, consistent with the usually quoted value for the 
incompressibility. 



B. Internal energy, specific heat and the ideal Bose gas 

The energy and specific heat are important quantities also, especially when the volume 
is fixed and the temperature is changing. The internal energy decomposes into low and high 
temperature components just as the pressure does, namely 

U = ^U„r^[l-^)ui,, (24) 

where U^r^ = Uur — Aeq (jf~^ "^{x+Thy simply the cutoff corrected low temperature 
energy. Figure |](a) plots this behavior for a number of densities, which reveals little except 
for the monotonic increasing nature of the energy. 
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The specific heat on the other hand is more reveahng. Using Cy = one arrives 

at 

Cv , ,d , , , , , 2kBT ( To 



kB " "1 ' ^ £o VT + Tc 

and Fig. |^(b) illustrates the behavior of Cy for d = 3, = k^^"^- The ideal gas limit is seen 
in the first term The low temperature Fermi gas result is contained in the second 

term {A — {m)){2kBT /eq). The last term involving the multiplicity fiuctuations gives rise 
to a peak in Cy/kB when d > 2. 

The variance induced peak is quite interesting. As we shall see in section p!V A|, it is 



analogous to condensation in an ideal Bose-Einstein gas, a critical transition which has been 
extensively studied. 

C. The entropy of nuclei 

The entropy of fragmenting nuclei can be obtained from the relation S = ^kBT In Z a 
which gives 

When (m) ^ 1, the entropy is that of a nearly degenerate Fermi gas of A nucleons confined 
to a cluster of size A, plus a contribution from the motion of this cluster which is coupled 
to a heat bath in the canonical ensemble. The nearly degenerate Fermi gas entropy is given 
by Eq. (^) where cutoff effects are excluded. Including it gets 

kB 2 \ eo I \ 2{T + T„)) ^ ' 

The cutoff correction is of the same order as the next higher order correction. When (m) ~ A, 
S reduces to the Sackur- Tetrode law of Eq. In the intermediate regime, their is a critical 
point where a latent heat must be overcome, causing a change in volume. 
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Figure ^ is a plot of S/ksA vs. ksT. In the plot, the contribution from the thermal 
motion of the largest cluster has been subtracted out so that S —>■ as (m) —>■ 1. This is 
equivalent to requiring the total momentum of the system to be zero, a correction consistently 
made throughout these calculations. 



IV. CRITICAL POINT BEHAVIOR AND PARALLELS WITH OTHER MODELS 

The thermodynamic functions strongly suggest a critical point, which we determine 
explicitly in this section. To start, let us work out some properties of the cluster dis- 
tributions in the grand canonical ensemble. Here the partition function is given by 
Z = J2A^'^ZAix,y) = ex'pJ2ki^/y)^''/k'^ where Pk = k'^ and z = z'y = e^^'^^'^ is the 
chemical potential. In this limit (n^) is simply 

k 

ink) = . (27) 
yk^ 

Since A = J2k f^i^k), we have 

^ = J2^'k'-^ = 9r-iiz), (28) 

^ k=l 

with x' = x/y, x,y as defined earlier and gn{z) = Z]fc>o ^""'-2^- Notice that gn^i{z) = 
zdgn/dz and that for 2; ^ 1, 

gniz) = u-'T{l - n) + E C(n - ^)(-l)'^ , (29) 

k>o 

where z/ = — In result due to London [p^ 



The LHS of Eq. (pSj) is always finite in the thermodynamic limit. If z > 1, the RHS 
always diverges in the infinite limit A ^ 00, and the equation cannot hold. At 2; = 1, 
J2kk^'^ = ('(''" the Riemann zeta function, which is finite only if r > 2. UA/x' < C(r— 1) 
then a z < 1 can be found that satisfies the last constraint equation including all terms in 
the sum. Otherwise, the grand canonical expression breaks down, which defines a critical 
point 
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4 = C(r - 1) . (30) 



The failure of the grand canonical ensemble is due to the appearance of an infinite cluster, 
i.e. of size a A where a is nonzero even as A tends to infinity. More specifically, when A/x' 
exceeds ({t — 1), an infinite cluster exists, absorbing enough mass so that the remaining 
mass is distributed grand canonically with z = 1. If A/x' is less that C(r — 1), an infinite 
cluster does not exist. 

As we have already seen, the multiplicity and its fluctuations play an important role in 
the thermodynamic functions. So how does the critical point affect the multiplicity and its 
fluctuations? The multiplicity (m) = Y.ki'^k) is given simply by 



(m) 
A 



x' < x' 

(31) 

'^^^ ' X > X\. 



The variance of the multiplicity Var(m) = x^{m) is given by 



Var(m) 



x' < x' 

^ ' (32) 



A I x'gr{z) _ gr-i(z) I ^ / 

Figure |^(a),(b) plots Var(m) = (m^) — (m)^. vs. x' and (m) respectively for both the 
canonical and grand canonical solutions. The discontinuity in the slope of Var(m) at the 
critical point is apparent from the figure. The behavior of Var(m) vs. x shows a cusp-like 
behavior. For finite systems, Var(m) has a rounded peak. As we shall see in the next section, 
this behavior is related to the ideal Bose-Einstein gas condensation. 

Is this critical behavior consistent with the experimental situation in nuclear collisions? 
At the critical point x = Xc, where Var(m) has a cusp, (m) = (m)^ with 

{m)c _ C{r) . . 

The EOS collaboration experiment determined that for gold multifragmentation, the charge 
multiplicity (m)c = 26 ± 1. Using the above expression and Z = 79 for A, results in a 
critical exponent r = 2.262 ± 0.013. Using a percolation theoretic analysis, they arrived at 
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a somewhat different r = 2.14 ± 0.06. The two results are sufficiently close to suggest that 
the connection between the critical multiplicity and critical exponent r |^0| may not differ 
greatly from this simple model. 



A. Ideal Bose-Einstein gas laws 

Let us now compare these results with the ideal Bose gas of A particles. The results 
can be found in Huang for the case d = 3 dimensions, and generalizing to arbitrary 
dimensions is fairly straightforward. 

In the grand canonical ensemble, the ideal Bose gas of A particles moving in a volume 

V at a temperature T has a fugacity z = exp{fi/kBT} determined by 

A 1 _ 1 

The critical temperature for a particular density and temperature occurs when 2; — > 1, 

V —>■ 00, which implies 

A 1 

Pc=y = JdddM^) (35) 

since above the critical point, the second term can be neglected in the infinite volume limit. 
Recalling that x = V/X^, this condition is identical to the fragmentation case if we replace 
X by x' and r by 1 + d/2. 

The energy, specific heat, pressure and incompressibility above and below the critical 
point are given by 

{ d{d+2) V / N d? 9d/2i^) r T 

Cv ^ I ^^xs9i+d/2{z} - tJ^JTA^ i > ^ ^^^^ 
ks 1 d{d+2) V ^ /I A T ^ T ' 

[ 4 Ad5'l+rf/2UJ 1 < J-c 

^^^¥..W^)T>T. ^^^^ 
^(?i+rf/2(l) T < T, 
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I T < T, 

Both above and below the critical point these results agree with the fragmentation case 



y = 1, r = 1 + d/2 if we apply Eqs. and (H). 



This agreement can be understood by noting that the weight given to the ideal Bose gas 
in d dimensions in a Mayer cluster expansion is given by Eq. (|15D with parameter vector 



^fc = + \ (40) 

At large x (i.e above the critical point) the second term is negligible and the model reverts 
to the one considered here. Below the critical point, the large cluster formed contributes 
little to the thermodynamic functions, allowing agreement to continue. 

The nature of the transition is first order as expected for a Bose-Einstein like condensa- 
tion. The particles which have accumulated into the zero momentum mode contribute no 
pressure and have no volume or entropy in the thermodynamic limit. The pressure due to 
the balance of particles is constant for a given temperature (i.e. independent of volume). A 
latent heat per particle of 

accompanies the transition and of course is the heat T/\S which must be spent in moving 
a particle from the vapor phase which has entropy to the condensed phase which has no 
entropy. The latent heat of Eq. (^l]) is the result of the change in volume along the constant 
pressure Maxwell line at a fixed T. The volume change is 



A 1 



= Ti—^, ■ (42) 



The entropy change can be obtained from the Clausius-Clapeyron equation dPc/dT = Sc/Vc 
to give 

Sc _ d + 2 ^ ^i+d/2(l) . s 

^ " ^^l^M^ ' ^^^^ 
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which is consistent with L = ksTSc as expected. The fraction of the number of particles in 
the Bose condensate (zero momentum mode) is 

for T < Tc. 

In nuclear fragmentation, this condensation is also present, but the clusterization occurs 
in real space instead of momentum space. Rather than the number of particles in the zero 
momentum mode signaling a phase transition, the number of particles in the largest cluster 
plays the same role. 

B. Critical point behavior and percolation theory 



Percolation models |^ have been used to describe nuclear fragmentation with percola- 



tion clusters corresponding to nuclear clusters arising from the collision process. The perco- 
lation cluster distribution is given in terms of the percolation probability p and percolation 
threshold probability above which an infinite cluster exists. Specifically, 

(rik) = ^f{ip-Pc)/Pc) , (45) 

where f{x) is a scaling function. Percolation models have also suggested useful methods of 
analyzing nuclear fragmentation data such as Campi plots 



The model in section has some features that are similar with the percolation model. 
First, only one parameter x' = x/y describes the distribution, the analog of p in percolation 
theory. Second, below a critical point x'^, the fraction of particles found in the largest cluster 
is finite, even if the number of particles tends to infinity, a situation familiar from percolation 
theory where for p > p^ a.n infinite cluster exists, while for p < p^ it can not exist. For the 
model of section |T| this can be seen by analogy with the Bose condensate, where a finite 
fraction of particles end up in the zero momentum state. In this nuclear fragmentation 
model, a finite fraction of the mass ends up in the largest cluster. Specifically if fcmax is the 
size of the largest cluster in each event, then in the infinite A limit 
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(^max) 

A 



(m) . 




(46) 



This threshold behavior of (/cmax)/^ as expected exactly parallels Bose-Einstein conden- 
sation. Figure |^ shows the behavior of (A^max) and its variance for some a range of x' for 
Pk = k^^"^- Due to scaling behavior (/^max) for various A's approach the single curve shown 
above. 

At the critical point, (n^) ~ A;"^ where r > 2 can be chosen when specifying the original 
weight. The critical exponent r for percolation is 2.21. 



V. CONCLUSION AND SUMMARY 

In this paper, hadronic matter under some extreme conditions was investigated using an 
exactly solvable model which can describe both the fragmentation of nuclei and the thermal 
properties in a unified way. Three distinct types of behavior characterize hadronic matter. 
At low temperatures, the fermionic aspects of the system are apparent through such quanti- 
ties as Fermi degeneracy pressure and Fermi energy. At high temperature and low density, 
ideal gas laws are shown to characterize the system. Between these two limits, a regime 
exists in which the behavior of the fragmenting matter resembles an ideal Bose gas with its 
associated Bose-Einstein condensation. The consequences of Bose-Einstein condensation on 
the properties of hadronic matter were studied. Simple expressions were obtained for the 
nuclear equation of state, heat capacity, compressibility, entropy, energy and fragmentation 
yields which encompass these three regimes. The importance of multiplicity fluctuations in 
statistical models of nuclear fragmentation and its effect on various thermodynamic prop- 
erties is stressed. Some recent data from the EOS collaboration on critical point behavior 
of hadronic matter is analyzed and simple expressions from the exactly solvable model are 
obtained which can account for some features of the data. 
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APPENDIX A: MULTIPLICITY EXPECTATION VALUES 



Given a function /(m) where m is the multiplicity, it is often useful to compute 
x-^{f(m{x,y))). If /(m) is a polynomial, the following identities can be used to compute 

d 

x—{[m]k) = {m[m]k) - {rn){[rn]k) , 
d 

x—{m^) = {m''^^) - {m){m^) , 

kJJu 

d 

x—{m)k = {m)k+i - k{m)2{m)k-i ■ 

where [x]n — x(x-l)....(x-n-|-l) and (m)„ = nth cumulant moment of m. The first is easily 
derived from the fact that ([m]fc) = '^^^ second identity can be derived from the 

first by expanding (m*^) = Y,j=o where |^| is the Stirhng number of the second 



is the 



kind, combined with the additional identity ([m]^) = Yl'i=o(—^)'' ^ ^- {n^'^), where 
unsigned Stirling number of the first kind, and the orthogonality relations between the two 
Stirling numbers. The third can be derived from the second using the binomial expansion 
of (m — {m)Y. 

These three identities are supplemented by three additional identities with a;^ — > —y-§-- 
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FIGURES 

FIG. 1. Pressure (a) and (in) compressibility (b) of the model at various temperatures. 

FIG. 2. Gibbs free energy per particle of the model as a function of pressure at various tem- 
peratures. 

FIG. 3. Internal energy (a) and specific heat (b) per particle of the model at various densities. 

FIG. 4. Entropy per particle of the model at various densities. 

FIG. 5. The variance in the nuclear multiplicity of the model as a function of x' /x'^ (a) and 
{m)/A (b). 

FIG. 6. Size of the largest cluster (a) and its variance (b) of the model as a function of x' /x'^. 
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